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We present a novel experimental scheme enabhng investigation of transient exotic spin couplings. 
The scheme is based on synchronous measurements of optical-magnetometer signals of several de- 
vices operating in magnetically shielded environments in distant locations (> 100 km). Although 
signatures of such exotic couplings may be present in a signal of the single magnetometer, it would 
be challenging to extract them from noise. With correlation measurements of signals from the mag- 
netometers, not only the effects can be identified but their nature may also be investigated. The 
ability of the network to investigate physics beyond the Standard Model is discussed by considering 
the spin coupling to stable topological defects (e.g. domain walls) of axion-like fields. It is shown 
that the network consisting of sensitive optical magnetometers is capable to probe an axion-like-field 
parameter space unconstrained by other experiments. 



I. INTRODUCTION 

Among all magnetometric techniques, optical mag- 
netometry P, [2] presently offers the possibility of the 
most sensitive magnetic-field measurements [3|. In- 
trinsic sensitivity of optical magnetometers (OMAGs) 
to spin dynamics does not limit their application to 
magnetic-field detection but also enables investigation 
of other spin interactions, including non-electromagnetic 
ones (see Ref. and references therein). In particu- 
lar, OMAGs can be applied to probe couplings between 
spins and hypothetical fields not predicted by the Stan- 
dard Model. Such exotic fields are postulated by a vari- 
ety of theories [l-[l^ . One manner in which they could 
manifest themselves on Earth is as transient events. A 
particular example would be transient coupling of spins 
to certain constituents of dark matter (DM) and dark 
energy (DE) 

Most experimental DM searches aim at direct detec- 
tion of some variety of particles that feebly interact with 
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ordinary baryonic matter, e.g.. Weakly Interacting Mas- 
sive Particles (WIMPs) or axions [31 • Until now, how- 
ever, all the searches have produced only upper limits on 
the interaction strength between DM and ordinary bary- 
onic matter. Over the years, however, alternative candi- 
dates for DM have been proposed. For example, if DM 
consists of light axions or axion-like particles, it behaves 
more like a coherent field than a collection of uncorrelated 
particles [Hill]. In some theoretical scenarios, because 
the vacuum energy of the axion field is non-zero, the 
field oscillates at a specific frequency and hence it would 
not produce static effects on baryonic matter. Such DM 
scenarios also tend to be generated by stable topological 
defects ll 4[l7H 23j. e.g., an axion-like field with a domain 
structure fl2| . When the Earth crosses one of the domain 
walls (DWs) separating regions with different vacuum 
expectation values of the axion-like field, a torc^ue can 
be exerted on leptonic or baryonic spins. Such a DW- 
crossing event could lead to transient sig nal detectable 
with modern state-of-the-art OMAGs [l^l- Based on as- 
tronomical constraints, however, one can show that wall- 
crossing events are rare and brief [l2j |. so the major issue 
becomes separation of the transient signals induced by 
the DW crossing from transient signals generated by en- 
vironmental and technical noise. Reliable rejection of 
OMAG's transient signals due to other effects requires 
development of a new approach. 

In this paper, we demonstrate the principles of a new 
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technique for detecting transient signals of exotic ori- 
gin using a global network of synchronized OMAGs (the 
Global Network of Optical Magnetometers for Exotic 
physics, GNOME). Although the network may be used 
for detection of a variety of transient interactions herald- 
ing physics beyond the Standard Model, here, for con- 
creteness, we focus our considerations on the transient 
effects induced by crossing through the DWs of axion-like 
field. We demonstrate that application of the GNOME 
enables probing presently unconstrain parameters of the 
field. 

The article is organized as follows. First, a general 
discussion of OMAGs with a special attention drawn to 
their most important characteristics relevant for detec- 
tion of exotic transient events is provided. Next, we re- 
view the theory of domain walls of axion-like fields rele- 
vant to their detection by the GNOME (Sec. Hill- Section 
IIVI discusses a first demonstration experiment using syn- 
chronously detected signals of two OMAGs separated by 
about 9000 km. These two magnetometers form the first 
elements of the envisioned future GNOME. The princi- 
ples that form the basis of the GNOME data analysis are 
outlined in Sec.|Vl Finally, the prospects of the GNOME 
are discussed in Sec. IVII and conclusions are drawn in 
Sec.lvni 



II. CHARACTERISTICS OF OMAGS 
RELEVANT FOR THE DETECTION OF 
TRANSIENT EFFECTS 

Detection of transient events that weakly perturb 
atomic spins requires OMAGs with specific characteris- 
tics. In particular, a suitable device needs to have high 
enough sensitivity to detect small changes in spin dy- 
namics. Moreover, its response to abrupt changes of spin 
behavior needs to be suitably fast not to distort or av- 
erage out the signals, so that their time-domain signa- 
ture can be reliably understood and compared between 
different GNOME sites. High sensitivity and high band- 
width, however, may not always be compatible with each 
other. Typically, high-sensitivity OMAGs also have char- 
acteristically slow spin relaxation times, which typically 
leads to narrower band widths, whereas high bandwidth 
devices typically have faster spin relaxation times, which 
degrades sensitivity. Below, we discuss fundamental and 
technical limitations of OMAG sensitivity as well as fac- 
tors determining the bandwidth of the devices. Discus- 
sion of other characteristics of the magnetometers, rele- 
vant for detection of transient signals, is also provided. 

In OMAGs, detection of magnetic fields occurs in a 
three-phase process. First, atoms are optically pumped; 
next, they evolve under the influence of external fields; 
finally, their quantum state is detected with light [63j . 
While this scheme allows for the most sensitive measure- 
ments of magnetic fields, it also sets a fundamental limit 
on the sensitivity of OMAGs. The limit results from the 
quantum nature of photons and atoms and the coupling 



between them. In the most general form, the fundamen- 
tal limit of the sensitivity SBf = ^ SBl^ + SB'^f^ + 6B^^, 

where 6Bat is the magnetic-sensitivity limit due to spin- 
projection noise (SPN), SBph is the limit related to pho- 
ton shot noise, and SB^a is the limit associated with back- 
action of the probe light on the atoms. The projection 
noise originates from the Heisenberg uncertainty princi- 
pie SF^SFf > |([F„F,])|V4 = ?i'(J^fc)V4, where F.^-fc 
are three components of the spin F and [,] denotes the 
commutator. When this relation become an equality, the 
SPN-limited magnetic-field sensitivity SBat may be writ- 
ten as [2|] 

giJ-B V NatT2T 

where Nat is the total number of atoms involved in 
the light-atom interaction, T2 is the transverse spin- 
relaxation time, T is the duration of the measurement, g 
is the Lande factor, is the Bohr magneton, and we use 
natural units {Ti — c = 1) [62|. Equation ([T]) reveals two 
strategies to improve the sensitivity of OMAGs. The first 
consists in prolonging the transverse spin-relaxation time 
T2, for example, by containing the atoms in a glass cell 
with antirelaxation coating that preserves spin polariza- 
tion upon atomic collisions with cell walls or introducing 
a buffer gas with a low spin-depolarization cross-section 
into the cell to limit diffusion to the walls. The second 
approach relies on increasing the number of atoms Nat- 
Both approaches are used in OMAGs and in fact, herein 
we discuss experimental results obtained with OMAGs 
exploiting both methods fSec. |V|. 

The second contribution 5Bph to the fundamental mag- 
netometric sensitivity 5Bf is associated with the fluctu- 
ation of the number of photons in the light beam used 
for probing the spins. Due to Poissonian statistics of 
photons, the intensity and polarization-state of light can 
only be determined with a finite precision cx (TVp/ir)"^/^, 
where Nph is a number of photons of the probe beam 
hitting detector per unit time. This sets a limit on the 
precision with which the spin state can be determined 
and hence the limit on the magnetometric sensitivity. It 
is important to note, however, that this contribution can 
be reduced by detuning the probe light frequency away 
from a resonant optical transition and simultaneously in- 
creasing the light intensity. The photon shot-noise lim- 
ited sensitivity improves due to the increase in Nph, and 
the probe light only weakly affects the medium while the 
state can still be efficiently determined (absorption on an 
isolated transition scales as I/A^ and dispersion as 1/A, 
where A is detuning). Hence, the contribution of the 
photon shot noise to the total magnetometric sensitivity 
5B} may be reduced so that SBph ^ SBat- 

The last source of fundamental noise comes from the 
Stark shift of energy levels induced by quantum fluctu- 
ations of light intensity and polarization (back action) 
[2^ ; fluctuations of energies of magnetic sublevels intro- 
duce uncertainty in spin precession and hence limit mag- 
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netometric sensitivity. Yet, there are means of reducing 
or eliminating backaction [23, [13, HI] . One approach is 
based on the same principle as used to reduce 5Bph, i.e., 
detuning the probe- light frequency far from any resonant 
transition [29]. This approach allows reduction of SB^a so 
that, under optimized conditions, Eq. ([T]) accurately rep- 
resents the fundamental sensitivity limit of the magnetic- 
field measurements with OMAGs {5Bf w SBat)- 

Typically, on top of fundamental noise, there is techni- 
cal noise SBf. For example, optical noise may be induced 
by mechanical vibration of optical elements or air turbu- 
lence in the probe-beam optical path. Electronics used 
in light detection can also contribute to technical noise. 
With appropriate experimental measures, the influence 
of the noise on overall magnetometric sensitivity may be 
reduced but it cannot be completely eliminated. Thus, 
in many cases, it is a significant (sometimes dominant) 
contribution to OMAG sensitivity. 

A different source of noise originates from uncontrol- 
lable magnetic fields. Such fields result in random re- 
sponse of OMAGs thus reduces the sensitivity of an 
OMAG to non-magnetic interactions affecting the atomic 
spins. In the case of OMAGs enclosed inside a mag- 
netic shield, a common configuration for precision mea- 
surements, uncertainty in spin-dynamics measurements 
may be introduced by external magnetic fields penetrat- 
ing into the shield [65|, thermal currents induced in the 
shield material, and instability of the current source used 
for generating magnetic fields. The noise may be reduced 
by application of active cancelation of the field outside 
the shield and/or used of low electric conductivity high 
magnetic susceptibility shielding materials [g^. Another 
common approach is to employ comagnetometry tech- 
niques, where the magnetic field is measured by multiple 
species expected to have different couplings to the exotic 
fields, allowing subtraction or cancelation of magnetic 
field noise. 

Although OMAGs do not have intrinsic 1/f noise, ex- 
istence of technical noise suggests an advantage of de- 
tection of optical signals at frequencies higher than 1//- 
"knee" . This may be achieved either by modulation of 
the probe light, i.e., by application of intensity, frequency, 
or polarization modulation, and phase-sensitive detection 
of the signal, or by operation of the devices in non-zero 
magnetic fields B h/T^gnB- In the latter case, the 
output signal of the magnetometer is modulated at the 
Larmor frequency lol = gf^sB/h or a multiple thereof, 
which enables filtering of the low-frequency noise. To 
detect such higher-frequency signals, however, OMAGs 
with broad dynamic ranges are required. 

Operation at non-zero magnetic fields raises another 
important issue in magnetic-field detection. Optical mag- 
netometers enable either scalar measurements, where the 
device response depends on the magnitude of a magnetic 
field, or vector measurements, where it is determined 
by specific vector components of the field. However, 
even scalar magnetometers operating at non-zero mag- 
netic fields become primarily sensitive to the field changes 



along the dominant component of the field; transverse 
components of the field add as second-order corrections 
to the total-field magnitude B. Moreover, modulation of 
the magnetic field in three spatial directions enables a 
scalar magnetometer to detect the three vector compo- 
nents of the field [sO]. There also exist techniques en- 
abling conversion of a scalar magnetometer into a vec- 
tor magnetometer without the necessity of applying a 
modulated magnetic field [Slj. The ability to determine 
not only a magnitude but also the direction of the spin- 
coupling field may have implications for the envisioned 
detection of transient effects due to exotic interactions. 

Another characteristic of OMAGs, particularly impor- 
tant in detection of transient signals, is bandwidth. For 
typical OMAGs (see Ref. [3] and references therein), 
the response of the magnetometer to small field changes 
is equivalent to a response of a first-order low-pass fil- 
ter with the time constant T2 [s^]. Hence the natural 
bandwidth of such OMAGs is given by {2ttT2)~^, which 
for shorter measurement times, i.e., t < T2, takes the 
form (27rr)~^. OMAG bandwidth can be broadened 
by shortening T2, which can be, for example, accom- 
plished by increasing intensity of the probe light (power 
broadening). That increase of the magnetometer band- 
width often occurs at the cost of its sensitivity [Eq. ([T])] . 
Therefore, optimized operation of OMAGs requires a 
compromise between the two quantities. It should be 
noted, however, that application of quantum nondemoli- 
tion measurements enables to achieve sensitive magnetic- 
field measurements at high bandwidth [33, 34]. 

In order to detect transient spin couplings, the signal 
characteristics must fall into the detection capabilities 
of the OMAGs used. Table U summarizes characteris- 
tics of various OMAGs with potential applicability to the 
GNOME. Although they differ in various aspects, many 
of them have the potential to be successfully employed 
as GNOME sensors. 



III. THEORETICAL BACKGROUND 

A specific example of exotic spin coupling that may be 
detectable with the GNOME is the transit of the Earth 
through a domain wall (DW) of a light pseudoscalar 
(axion-like) field [l^]- Stable domain structure of axion- 
like fields is a consequence of certain Standard Model 
extensions (45l - l48l |. Domains form out of the initially 
random distribution of the vacuum expectation values of 
the axion-like field as the Universe expands and cools. In 
this scenario, DWs separate regions of space with differ- 
ent energy vacua [4§] . Importantly, based on astrophysi- 
cal constraints, only light axion-like fields can build DWs 
that persist to the present epoch [67] . 

A detailed theoretical background of the optical detec- 
tion of wall crossings is presented elsewhere [l^]. Here 
we only briefly review the concept. We start with con- 
sidering a hypothetical pseudoscalar field a(r) that per- 
meates the Universe and forms a domain structure. As 
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TABLE I: Various OMAGs characteristics important for detecting transient signals due to exotic spin couplings. SEf and SEd 
are the fundamental and demonstrated OMAG sensitivities in energy units. The names of the magnetometers indicate the type 
of the device. HFP is an abbreviation of hexafluorobenzene. 



Name 


Element (s) / 
Compound(s) 


SEf 
(^eV/\/Hz^ 


oE 


T2 (ms) 


Spin 
coupling 


Ref. 


SERF 


K+Na+^He 


2.9 X 10-^2 


3.1 X 10"2i 


100 


Total 


f3] 


/i-SERF 


Cs 


3.9 X 10"^° 


5.8 X 10^2° 


5 


Total 


[35] 


NMR-SERF hybrid 


pentane-HFB 


^Q-23 


3.5 X 10"" 


20000 


Nuclear 


[36] 


NMOR 


Rb 


5.8 X 10"^^ 


5.8 X 10"" 


300 


Total 


[37] 


FM/AM NMOR 


Rb 


4 X 10~^° 


1.2 X 10"^* 


300 


Total 


[38-40] 




K, Cs 


1.4 X 10"" 


5.8 X 10"" 


300 


Total 


[41,42] 




Cs 


7.2 X 10"" 


1.4 X 10"^* 


0.06 


Total 


[43] 


Hg EDM 


Hg 


2x10"^'' 


2.5 X 10^2^ 


100000 


Nuclear 


m 



shown in [l^, a specific realization of the field exist- 
ing between neighboring domains with different energy- 
degenerate vacua (with the DW centered at z=0) 



a{z) — 4ao arctan[exp(maz)], 



(2) 



where ao is the characteristic amplitude of the field and 
nia is the pseudoscalar-particle mass. Coupling between 
the axion-like-field gradient Va and the spin F aris- 
ing during the domain-wall crossing is described by the 
Hamiltonian 163 



F • Vc 



of 



(3) 



where /of is the effective decay constant in units of en- 
ergy, /ef depends on the atomic structure of the particles 
used in a specific OMAG and is a combination of elec- 
tron /c, proton /p, and nucleon decay constants f^. By 
substituting Eq. ([2]) into Eq. ([3|), the Hamiltonian Huw 
can be expressed using the field parameters iria and Uq 



2 aouia cos ip 
fci cosh{maz) ' 



(4) 



where ip is the angle between the spin F and the field 
gradient Va. 

Until now, only lower bounds on the electron /e, neu- 
tron /„, and electron fp proton decay constants have 
been established by astronomical observations (|/c,n,p| > 
10^ GeV). Below we show that the coupling can be fur- 
ther investigated using the GNOME. 

The thickness of the DW d is determined by the 
pseudoscalar-particle mass rua via 



(5) 



Consequently, the mass also limits the duration of the 
transient signal At = d/v±, where v± is the relative 
speed between the DW and the OMAG. At distances 
much larger than d, the walls may be characterized by 
the tension cr, which is the mass/energy per unit area of 



the DW. In the considered case, it can be written as a 
function of the field parameters 



(6) 



Tlie tension can be related to the DW energy den- 
sity pDW via pdw ~ '^/L, where L is the charac- 
teristic size of the domain. Importantly, the density 
Pdw needs to be smaller than the DM density pdm 
(pDM ~ 0.4 GeV/cm'^) or the DE energy density poE 
(PDE » 0.4 X 10"^ GeV/cm^). Determination of the 
tension also requires knowledge about the characteristic 
size of the domain L. Since it is not possible to deter- 
mine L without further assumptions about the specific 
mechanism of domain-structure formation, here we treat 
L as a free parameter and constrain it from an experi- 
mental perspective, i.e., the experimental feasibility im- 
plies that the average time T between two wall crossings 
should not be longer than 10 years. By taking into ac- 
count the speed of the solar system relative to the Galac- 
tic frame {v « 10^'^c), a DW-crossing event will occur 
within a time-span of 10 years if the domain size is less 
than 10"^ ly 

Combining Eqs. ([4]) and one obtains the formula 
for the effective decay constant /ot at the center of the 
wall (z = 0) 



/cf = cos p, 



(7) 



where SE is the sensitivity of an OMAG in energy units. 

Figure [T] presents the parameter space that can be 
probed with OMAGs. The shaded region indicates the 
parameter range that can be realistically probed with 
the OMAG of a sensitivity of w 3 x lO^^o qY/VBz and a 
characteristic DW-size of 10"^ ly. The plot indicates that 
detection of axion-like field domain structure with the 
DM energy density pdm is possible with the envisioned 
GNOME. Moreover, with higher sensitivity magnetome- 
ters (see Table |T|, the parameter space can be further 
explored enabling probing DWs with the density pdw, 
given by the DE density pde 
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FIG. 1: Parameter space of the axion-like field with a do- 
main structure that can be probed with the GNOME (green 
shaded region). The blue line limits the space due to the 
sensitivity of applied magnetometers (assumed here to be 
2.9 X 10~^° eV/\/Hz, which corresponds to a sensitivity of 
■j^Q-15 t/VHz of the cesium NMOR magnetometer). The two 
green lines arise from exiting limits from the astronomical 
observations [SQ] (the horizontal line at lO'"* GeV) and the du- 
ration of the experiment being less than 10 years (vertical line 
at 10~* eV). For the plot the DW energy density was assumed 
to equal the DM energy density (pdm ~ 0.4 GeV/cm^). 



Being able to detect a single lO-ms duration event in 
a 10-year time-span requires continuous and robust op- 
eration of the magnetometer on a comparable time scale. 
Figure[2]presents the average time between DW crossings 
as a function of the decay constant /of an OMAG with a 
sensitivity equal to 2.9 x 10^° eV/VSz (1 FT/Vb^)- The 




FIG. 2: Average time between successive DW crossings T as 
a function of the effective decay constant /ef (green shaded). 
The blue line is determined by a magnetometer sensitivity of 
2.9 X 10^° eV/\/HS corresponding to 10 T/^HS), while the 
meaning of green lines is the same as in Fig. [T] 



FIG. 3; The concept of the synchronized-magnetometer ar- 
rangement. OMAGs located at globally separated locations 
record signals with a time synchronization provided by the 
GPS. By synchronously detecting and correlating magne- 
tometer signals, transient events of global character may be 
identified. In particular, correlating signals of at least four 
devices enables detection of spatiotemporal character of the 
event. Here, two devices located in Krakow and Berkeley are 
used to search for transient signals induced by crossing of a 
DW of an axion-like field. 



In principle, the parameters of the model can be con- 
strained with a single magnetometer. A particular prob- 
lem for a search carried out with a single OMAG is the 
appearance of brief spikes in the OMAG signal related 
to technical noise or abrupt magnetic field changes. In 
a single device, rejection of these false-positive signals is 
difficult. At the same time, coincident measurements be- 
tween two or more instruments are helpful in rejecting 
such signals; they provide consistency checks, since a sig- 
nal would be expected to exist in all instruments whereas 
environmentally induced events are not typically corre- 
lated in the time window required for coincidence. Fur- 
thermore, information about a putative event such as its 
impinging direction can be determined by triangulation 
if several instruments (at least four) are taking data si- 
multaneously (see discussion in Sec. IVip. These features 
clearly show that synchronous operation of multiple syn- 
chronized, geographically separated OMAGs within the 
proposed global network may facilitate searches for such 
transient signals of astrophysical origin. 



IV. EXPERIMENTAL APPARATUS 



results show that the GNOME can probe significant re- 
gions of parameter space in as few as three days of contin- 
ues operation. The accessible parameter space expands 
rapidly with increasing duration of the measurements. 



The concept of the experimental apparatus is shown 
in Fig. [31 Since detailed description of the two OMAGs 
used in our experiment may be found elsewhere (ssl . l5l| , 
we only briefiy review the most important elements of 
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the devices. Both magnetometers use rubidium va- 
pors as the magneto-optically active medium. In the 
Krakow magnetometer, the atomic vapor is contained 
in a parafSn-coated evacuated cyhndrical glass cell with 
volume ~ 3 cm"^. The vapor cell is maintained at about 
50°C corresponding to an atomic density of roughly 
10"'^-'^ atoms/cm'^. The relaxation rate of the atomic 
ground state is 27r x 30 s~^, which yields a fundamen- 
tal sensitivity SBf of w 3 fT/VHz (spin-projection lim- 
ited) [6^. The second magnetometer (Berkeley) exploits 
a microfabricated vapor cell 5^ of a volume of 0.01 cm^ 
that is heated up to about 200°C. Operation in the 
spin-exchange relaxation free regime [53| allows elimina- 
tions of relaxation due to spin-exchange collisions, one 
of the main ground-state polarization-relaxation mech- 
anisms. Application of the technique allows one to ob- 
tain a ground-state relaxation rate of about 27r x 400 s~^, 
which in combination with 3-4 orders of magnitude higher 
density yields a similar sensitivity as for the other setup, 
i.e., about 1 fT/v^Hz (t^]. Both magnetometers thus have 
the capability to detect a DW crossing and probe the pa- 
rameter space. 

Both magnetometers are placed inside multilayer mag- 
netic shields made of /x-metal with the innermost layer 
made of ferrite [7l|. The shields reduce external mag- 
netic fields by a factor 10^. Inside the shield atoms are 
subjected to a stable, well-controlled magnetic field gen- 
erated by a set of three-dimensional magnetic-field coils. 
In the Krakow magnetometer a field with a magnitude 
of 10^'' T is applied, while at Berkeley the applied-field 
magnitude is ~ 5 x 10~* T. 

The outputs of the magnetometers are acquired us- 
ing custom-made devices based on Trimble Resolution-T 
GPS (Global Positioning System) time receivers. The 
data acquisition devices provides time markers separated 
by one second with a precision of about 80 ns synchro- 
nized with a quartz clock built into the devices. The 
acquisition devices can record simultaneously signals in 
four channels with at a rate of 1000 samples/s. Each 
one-second-long record is stored on a memory card with 
a header containing information on time, location, etc. 
Every record is transmitted to the computer (via serial 
port) where it is binned into groups of 10-1000 records 
(typically 2-minute long bins are generated). The data 
are stored with computers located at the respective loca- 
tions, and every 1-2 hours the information is exchanged 
between Krakow and Berkeley using File Transfer Pro- 
tocol (FTP). In this manner, the complete set of data is 
accessible in either of the locations. The detailed discus- 
sion of the GPS-synchronized data-acquisition devices is 
presented in Ref. j53 |. 



V. RESULTS AND DATA ANALYSIS 

Figure 3] presents magnetometer signals measured syn- 
chronously in the two locations (Berkeley, California, 
USA and Krakow, Poland) over a period of about 



0.30- 



0.25 



.ff 0.20- 



0.15- 



O.IOL 



Berkeley 



Kiakow 



0.9 



0.75 



0.6 5) 



0.45 



0.3 

1000 2000 3000 4000 5000 6000 

Time [s] 

FIG. 4: Synchronously detected signals of the two OMAGs 
located in Krakow and Berkeley. The signal in the Krakow 
magnetometer was measured at a field of 100 nT, whereas the 
Berkeley magnetometer operated at 50 nT. The data acquisi- 
tion rate was 256 samples per second. 



1.5 hours. Relative to the mean amplitude of the signal, 
the long-term drift of the Krakow magnetometer is larger 
than that of the Berkeley magnetometer. This is most 
likely attributable to instability of the laser frequency, 
which will be addressed in the future by implementation 
of laser-stabilization techniques. The Berkeley magne- 
tometer signal exhibits short duration (< 4 ms) spikes of 
relatively large amplitudes not observed in the Krakow 
signal. Auxiliary tests verified that the noise spikes orig- 
inated from electronic pick-up, a problem that will be 
addressed in the future. 

In many respects, identification of a DW-crossing event 
using the GNOME is similar to searches for gravitational- 
wave bursts with a system of long-baseline laser interfer- 
ometers [such as the Laser Interferometer Gravitational- 
Wave Observatory (LIGO), the Virgo detector, GEO 600, 
and TAMA 300] f55'|. Both types of experiments aim 
to identify and characterize transient signals and search 
for time-domain correlations between the transient sig- 
nals from different detectors. Importantly, the field of 
transient gravitational-wave astronomy has developed a 
variety of statistical methods to identify brief (duration 
< 1 s) signals correlated among different detectors but 
otherwise generic in noisy time series data. 

As a proof-of-principle demonstration, we have ap- 
plied one of the methods upon which such statistical 
analysis is based, the "excess power" statistic [s^, to 
the synchronous magnetometer data from the Krakow 
and Berkeley sites. The analysis is carried out as fol- 
lows. First, an estimation of the power spectral den- 
sity (PSD) over several continuous, overlapping segments 
of the data is made for the data from each individual 
detector. These spectra are updated at regular inter- 
vals and combined with previous measurements using a 
running median exponentially-weighted history (this is a 
filter that weights contributions to the moving average 
by a factor that decreases exponentially with the time 
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since the data was acquired). The PSD is then used to 
whiten the data, ideahy, producing a stream of Gaus- 
sian distributed, zero-mean, unity variance random vari- 
ables characterizing the data set. Whitening the data 
removes correlations between the different variables used 
to characterize the magnetometer signals. The stream of 
whitened data is then passed through a bank of band- 
limited filters producing several channels of filtered data. 
The filters are Hann windows in the frequency domain, 
which are themselves whitened with the PSD. The re- 
sulting streams are a discrete localization of the energy 
in the original data stream described by "tiles" having a 
time-frequency extent bounded by a bandwidth Ai// (the 
bandwidth of the filters) and duration 6t. The resulting 
tile is constructed so as to have iVDOF independent de- 
grees of freedom such that A^dof = 2Si'f6t [57^ . 
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FIG. 5: Time-frequency maps of outlying tile energies of the 
Berkeley (top) and Krakow (bottom) OMAG data. Each 
dot represents a tile characterized by its power-weighted cen- 
tral time (a;-axis), Fourier frequency (y-axis), and tile energy 
(color and dot size). The large markers indicate tiles with a 
normalized tile energy greater than 100. Large tile energies 
are likely caused by discontinuities (which cause short but 
broadband responses) in the data or environmental influences 
within the instruments themselves. 

The final product is a time/frequency map of tiles de- 
scribing the whitened signal energy (Fig. [SJ. Under the 
conditions of stationarity (the PSD does not fluctuate on 
the time scale of the estimation process) and Gaussian- 
ity (the data samples have a distribution matching the 
Gaussian distribution), the tile energies are distributed 
as a distribution with degrees of freedom as described 
above. Thus the significance of any tile's energy is well- 
understood and the statistical probability of outlier tiles 
can be measured. Correlated transient signals from dif- 
ferent magnetometers can be searched for by using the 
time / frequency tile maps to find overlapping events with 
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FIG. 6: Histogram showing the number of tiles with energy 
within a specific energy range for the magnetometer data de- 
picted in Figs. [4] and [5] 



related characteristics. Time shifts can be introduced 
into the data to check for correlated transient events with 
particular relative delays. 

Examination of the time/frequency map of whitened 
signal energy shown in Fig. [S] reveals that the Krakow 
instrument produces fewer "loud" transients (associated 
with high tile energies) than the Berkeley instrument, in- 
dicating the that Krakow data more closely follows the 
signal expected from Gaussian-distributed noise. The 
loud transients in both instruments are likely due to elec- 
tromagnetic environmental couplings or technical noise 
in the instruments themselves (laser light power and fre- 
quency noise, etc.). Further investigation is needed to 
identify the causes of the loud transients and find ap- 
propriate methods to eliminate as many sources of loud 
transients as possible. 

Figure [H] is a histogram of the tile energies from the 
trigger maps in Fig. [S] with the Berkeley magnetometer 
signal in blue and the Krakow magnetometer signal in 
green. One observes a roughly one order-of-magnitude 
excess of events in the Berkeley magnetometer signal 
compared to the Krakow signal for tile energies above 
^ 100. This indicates problems in the relative quality 
of the Berkeley data compared to the Krakow data, per- 
haps a result of environmental noise or technical noise 
(such as the electronic pick-up noise identified in auxil- 
iary measurements). Such histogram analysis is another 
useful way to characterize the performance of individual 
magnetometers in the GNOME. 



VI. PROSPECTS 

Comagnetometry, where the magnetic field is simul- 
taneously measured with multiple atomic species or de- 
vices, is a widely used technique in precision mea- 
surements searching for anomalous spin-dependent ef- 
fects. Comagnetometry with different atomic species of- 
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ten takes advantage of the fact that the relative couphng 
strengths of an exotic field to electrons and nuclei are 
generally different from the relative coupling strengths 
of electrons and nuclei to magnetic fields. A particular 
example of a comagnetometry scheme that will be in- 
vestigated for possible use as a GNOME sensor is SERF 
comagnetometer similar to that described in Ref . [58]] . In 
contrast to the devices used in the demonstration exper- 
iment described in the present work, an additional noble 
gas (helium) is introduced into a vapor cell. When the 
noble gas has non-zero nuclear spin, the alkali and noble 
gas spins become strongly coupled through spin-exchange 
collisions [H, 113] . This coupling can be represented as 
the effective magnetic field BeS experienced by one spin 
species due to the average magnetization M of the other, 
due to enhancement of the alkali valence electron density 
at the noble gas nucleus, 

Beff = AM, (8) 

where A is a parameter determined by the particular 
properties of the alkali-noble gas spin-exchange [6l[ . The 
applied field B is tuned so that it approximately cancels 
Boff experienced by the alkali atoms. The alkali atoms 
are then in an effective zero-field environment, and be- 
cause the noble gas magnetization M adiabatically fol- 
lows B, transverse components of B are automatically 
compensated by Beg' to first order. Such cancelation only 
occurs for interactions that couple to spins in proportion 
to their magnetic moments, leaving the SERF comag- 
netometer sensitive to anomalous spin couplings to elec- 
trons and nuclei [59| . 

The energy resolution of the latest generation of the 
SERF comagnetometer, employing Rb as the alkali atom 
and 2iNe as the noble gas, is IQ-^^ eV/v1iz [Hi. This 
new scheme uses hybrid optical pumping of Rb via spin- 
exchange collisions with low-density, optically pumped 
K and off-resonant direct optical probing of Rb spins. 
This approach allows full optimization of both optical 
pumping and probing. Because of the relatively small 
gyromagnetic ratio of ^^Ne, the Rb-K-^^Ne SERF co- 
magnetometer has an order of magnitude better energy 
resolution for the same level of magnetic-field sensitivity 
as compared to earlier SERF comagnetometers, and may 
offer advantages in bandwidth. In the future, we plan to 
develop and optimize the SERF-based comagnetometer 
for measurements of exotic transient effects. 

Independently from the development of SERF-based 
comagnetometer, the other magnetometer types will be 
developed as potential GNOME sensors. In particular, 
we envision using sensors that monitor evolution of var- 
ious types of spins (proton, neutron, electron). This 
would add another dimension to our investigations by 
studying infiuence of exotic coupling to various funda- 
mental particles. 

Another important work envisioned for a future exper- 
iment is correlation of the magnetometer readouts with 
environmental parameters (e.g., magnetic field outside 
the shield, temperature, etc.). This is motivated by the 



fact that despite magnetic shielding, there will inevitably 
be some level of transient signals and noise associated 
with the local environment (and possibly with global ef- 
fects like the solar wind, changes to the Earths magnetic 
field, etc.). The environmental-condition data will allow 
for exclusion/ vetoing of data with known systematic is- 
sues. 

Further step in reducing the infiuence of magnetic 
fields on the operation of GNOME is application of Su- 
perconducting Quantum Interference Device (SQUID) 
[6 21 magnetometers as sensors operating in addition to 
OMAGs inside the magnetic shields. While the SQUID 
magnetometers are characterized with magnetometric 
sensitivity comparable to that of OMAGs, they are not 
sensitive to exotic spin coupling. Thus, they can be used 
for vetoing false-positive transient signals. 

Ultimately, the GNOME will consist of at least five 
OMAGs. Four devices will be used for the detection of 
a DW and of its geometrical properties. Any additional 
magnetometer would increase the sensitivity of the net- 
work. An independent OMAG will serve as cross-check 
to verify if, based on predicated DW event, a transient 
signal arises in the magnetometer in a narrow temporal 
window. 



VII. CONCLUSIONS 

In this paper, we presented a new experimental scheme 
enabling investigations of transient exotic spin couplings. 
It is based on synchronous operation of globally separated 
optical magnetometers enclosed inside magnetic shields. 
Correlation of magnetometers' readouts enables filtering 
local signals induced by environmental and/or technical 
noise. Moreover, application of vetoing techniques, e.g., 
via correlation of optical-magnetometer readouts with 
signals detected with non-optical magnetic-field sensors, 
enables suppression of influence of global disturbances of 
magnetic origins, such as solar wind, fluctuation of the 
Earth's magnetic fleld, on the operation of the magne- 
tometers. In such an arrangement, the network becomes 
primarily sensitive to spin coupling of non-magnetic ori- 
gins, thus it may be used for searches of physics be- 
yond the Standard Model. A specific example of such 
searches was discussed here by considering coupling of 
atomic spins to domain walls of axion-like fields. It was 
demonstrated that with modern state-of-the-art optical 
magnetometers probing a significant region of currently 
unconstrained space of parameters of the fields is feasible. 
The preliminary results obtained based on synchronous 
operation of two magnetometers located in Krakow and 
Berkeley were presented and future plans for the network 
development were outlined. 
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[63] The phases may either be separated in time or may occur 
simultaneously. 

[64] Although for a minimum uncertainty state, the sensitiv- 
ity SBat scales as 1/ y/Nat, for a quantum system with en- 
tanglement the scaling could be stronger (SBat oc 1/Nat)- 
In principle, it should allow for a large improvement in 
magnetic-field sensitivity of OMAGs. Unfortunately, the 
entangled states are fragile and rarely improvement of 
the sensitivity SBat below Eq. ([T]) is observed [2^. [25l]. 

[65] Note that the sensitivity of the magnetometer is at a level 
of 10~^^ T/VHz or better, while the Earth magnetic field 
is ^ 4 X 10~^ T. Hence even with a shielding factor of 
> lO'' for DC fields, the magnetometer is still strongly 
sensitive to the fluctuation of the external flelds. 

[66] In order to limit magnetic-field noise due to Johnson cur- 
rents in the shield, modern magnetic shields have the 
innermost layer made of ferrite. Although such mate- 
rial has significantly lower magnetic susceptibility than 
permaloys, in particular, ^-metal, they have orders of 
magnitude larger resistivity, which suppresses thermal 
currents in the magnetic-shield layer placed closest to 
a vapor cell (see Ref. Q and references therein) . 

[67] Formation of the wall from QCD-axion field would lead to 
disastrous cosmological consequences due to the excessive 
energy stored in the walls. 

[68] Note that in general, the higher-order couplings to the 
spins may be considered (see Ref for further details). 

[69] Note that the demonstrated sensitivity of the magne- 
tometer SB is lower. 

[70] The demonstrated sensitivity of the device SB is roughly 
10 fT/VHi. 

[71] While magnetic shields do not screen exotic interactions, 
their role in searches for transient exotic couplings require 
more thorough investigations in the future. 



